Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



IMAGING FABRY-PEROT SPECTROSCOPY OF NGC 5775: KINEMATICS OF THE DIFFUSE IONIZED GAS 

HALO 

George H. Heald and Richard J. Rand 

University of New Mexico, Department of Physics and Astronomy, 800 Yale Boulevard NE, Albuquerque, NM 87131 

Robert A. Benjamin 

University of Wisconsin — Whitewater, Department of Physics, 800 West Main Street, Whitewater, WI 53190 

Joseph A. Collins 

University of Colorado, Center for Astrophysics and Space Astronomy, Campus Box 389, Boulder, CO 80309 

AND 

Joss Bland-Hawthorn 

Anglo-Australian Observatory, P. O. Box 296, Epping, NSW 2121, Australia 
Draft version February 5, 2008 

ABSTRACT 

We present imaging Fabry-Perot observations of Ha emission in the nearly edge-on spiral galaxy 
NGC 5775. We have derived a rotation curve and a radial density profile along the major axis by 
examining position-velocity (PV) diagrams from the Fabry-Perot data cube as well as a CO 2-1 data 
cube from the literature. PV diagrams constructed parallel to the major axis are used to examine 
changes in azimuthal velocity as a function of height above the midplane. The results of this analysis 
reveal the presence of a vertical gradient in azimuthal velocity. The magnitude of this gradient is 
approximately 1 km s -1 arcsec -1 , or about 8 km s _1 kpc -1 , though a higher value of the gradient 
may be appropriate in localized regions of the halo. The evidence for an azimuthal velocity gradient 
is much stronger for the approaching half of the galaxy, although earlier slit spectra are consistent 
with a gradient on both sides. There is evidence for an outward radial redistribution of gas in the 
halo. The form of the rotation curve may also change with height, but this is not certain. We 
compare these results with those of an entirely ballistic model of a disk-halo flow. The model predicts 
a vertical gradient in azimuthal velocity which is shallower than the observed gradient, indicating that 
an additional mechanism is required to further slow the rotation speeds in the halo. 
Subject headings: galaxies: halos — galaxies: individual (NGC 5775) — galaxies: kinematics and 
dynamics 



1. INTRODUCTION 

Gaseous thick disks are observed in the Milky Way and 
in some external edge-on spirals. In the Milky Way, the 
layer of vertically extended ionized emission is known as 
the Reynolds Layer or the warm ionized medium (WIM). 
This phase consists of gas with T ~ 10 4 K an d n R ~ 0.1 
cm" 3 , and has a scale height of about 1 kpc (Rcvnolds 
1993). 

In external galaxies, these diffuse ionized gas (DIG) 
layers are observed to have widely vary i ng morphologies 
(e.g.. iRand. Kulkarni. fc Hesterl 1199(1 iDettmarl 1199(1 
Hoones. Walterbos. fc R andl 119991: iMiller fc Veilleuxl 
2003). Not all edge-on spirals show detectable ex- 
traplanar (i.e., above the layer of H II regions) DIG 
(EDIG) emission. The prominence of a galaxy's EDIG 
layer is now known to correlate with tracers of star 
formation in the disk, such as the surface density of FIR 
emission, the star formation rate as deter mined by Ha 
luminosity, and the dus t tempera ture (e.g.. IRandlll996t 
iRossa fc Dettmarl 120001 l2003albD . Other gas phases, 
as well as dust absorption, have also been observed 
high above the midplane (see e.g.. ICollins et alJ 120001 
iLee et al.ll200lt lDettmarll2005l IFraternah et al.ll2005|) 

These observations are generally interpreted in terms 



of a star formation driven disk-halo flow such as t hat 
described by the chim ney lINorman fc Ikeuc h1l | 1989l) or 
galac tic fountain fe.g.. IShapiro fc Field! Il976t iBregmanl 
1980) model, and we adopt that viewpoint in this pa- 
per. Rising hot gas may originate in chimneys but warm, 
swept up ambient gas may also be pushed upwards to 
significant heights. The evolution of such gas and the in- 
teraction of different extraplanar phases is not at all well 
understood. However, the observations described above 
strongly suggest that such flows exist. 

A likely effect of disk-halo cycling is redistribution 
of gas in t he disk. The mass flux rates estimated 
thus far (cf. iBregman fc Pildisill99^: IWang et all 119951 
IFraternah et alJl20 02) imply that this process is respon- 
sible for moving a large amount of gas. Whether the cy- 
cling process is dominated by ballistic motion, (magneto- 
) hydrodynamic effects, or both, is still an open question. 
But the cycling may affect the distribution and rate of 
star formation in the dis k, and may result in e xtrapla- 
nar star formation (e.g., iTullmann et alJ 12003(1 . That 
a disk-halo flow is responsible for some of the observed 
high-velocity clouds (HVCs) is also a possibility (e.g., 
Wakk er"fc van Woerdenlfl997j) . For these reasons, it is 
important to understand the kinematics of this disk-halo 
cycling. 
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Studies of extraplanar HI emission have so far pro- 
vided the most spatially complete velocity information 
about gas which may be participating in disk-halo in- 
teractions. Halo gas is seen to lag the underlying disk 
rotat i on in NGC 891 l l Swaters. Sancisi. fc van der Hulst 
1997} iFraternali et all I2005T). NGC 5775 llLee et all 
2001), and NGC 2403 (Schaap, Sancisi. fc Swatersil200C 



TABLE 1 
Galaxy Parameters for NGC 
5775 



Fraternali et al-lfeOOlj) . This result is expected in a disk- 



halo flow model because as gas is lifted into the halo, it 
feels a weaker gravitational potential, migrates radially 
outward, and thus its rotation speed drops in order to 
conserve angular momentum. Since DIG has turned out 
to be an excellent tracer of disk-halo cycling, its kine- 
matics may shed significant light on the nature of such 
flows. 

Limited studi es of EDIG kinematics have been per- 
formed to date. IFraternali. Oosterloo. fc Sancisil l|2004D 
use spectra from slits along the major and minor axes 
of NGC 2403 to detect a rotational lag in the extra- 
planar ionized gas component. Spectra have also been 
obtained along slits perpendicul ar to the major axes of 
NGC 891 (one slit: IPlmlll997h and N GC 5775 (three 
slits: lRarT^l2(TfKl iTiillmann et al.ll200fl() . In both galax- 
ies, the mean velocities are seen to approach the sys- 
temic velocity as height above the midplane (z) increases. 
Models representing two distinct physical regimes have 
been generated to try to understand the apparent drop- 
off in rotation speed with z in NGC 891 and NGC 
5775. Purely hydrodynamical models are desc ribed 
and compared t o exis ting data by iBeniamml l)2000D and 
iBarnabe et al.l (120051) : a purely ballis ti c mo del is pre- 
sented bv iCollins. Benjamin, fc Randl (|2002, hereafter 
CBR). 

The previous studies of EDIG kinematics were limited 
by poor velocity resolution and the one-dimensional spa- 
tial coverage of long-slit spectroscopy. The mean veloci- 
ties used are affected by both the rotation speed and the 
distribution of emission along the line of sight. A full 
analysis of how gaseous halos rotate requires that these 
effects be separated as has been done with H I observa- 
tions of NGC 891 (e.g.. IFraternali et alJl2005l) . and this 
can only be achieved with high spectral resolution and 
two-dimensional spatial coverage. The work described in 
this paper represents an extension of the previous optical 
studies, in that high spectral resolution velocity informa- 
tion is obtained and analyzed for the full two-dimensional 
extent of NGC 5775. 

NGC 5775 is classified as an SBc galaxy. It is un- 
dergoing an interaction with its neighbor, NGC 5774, 
and a tidal s tream of HI connec ts the two galaxies (see 
Iirwinlll994fi . ICollins et all (j2000D cite a high far-infrared 
luminosity determined from the IRAS satellite, Lfir = 
7.9 x 10 43 erg s _1 , and a "far-infrared surface brightness" 
of L Fm /Z>| 5 = 8.4 x 10 40 erg s" 1 kpc~ 2 (where D 25 is 
the optical isophotal diameter at 25 th magnitude), whic h 
is typical of mild starbursts ijRossa fc Dettmar1 l2003a). 
The EDIG layer is bri ght, with a large sca le height and 
filamentary structures ijCollins et al.ll2f)00[) . A summary 
of galaxy parameters for NGC 5775 is presented in Table 
□ 

This paper is organized as follows. We describe the 
observations and the data reduction steps in In [0 
we present an analysis of disk and halo rotation. We 



Parameter 



Value 



RA (J2000.0) 
Decl. (J2000.0) 
Adopted Distance 51 
Inclination 
Position Angle 
Systemic Velocity 



14 h 53 m 57?57 
03 d 32 m 4Wl 
24.8 Mpc 
86° 
145.7° 

1681.1 km s" 1 



Not e. — All values are from llrwinl 
<199l . 

a Assuming Hq = 75 km s _1 Mpc - 1 



compare the data with the ballistic model of CBR in 21 
and conclude the paper in 33 

2. OBSERVATIONS AND DATA REDUCTION 

Data were obtained during the nights of 2001 
April 11-13 at the Anglo- Australian Telescope 
(AAT). The TAURUS-II Fabry-Perot interferome- 
ter, which is placed at the Cassegrain focus (f/8) of 
the AAT, was used in conjunction with the MIT/LL 
2kx4k CCD. Design, theory, and data reduction 
techniques of Fabry-Perot interferom eters are well 
described in the literature; see e.g. iBland fc Tullvl 
llT989Tl: [Ws. Shonhell. fc Bland-Hawthornl }2002) ; 
iGordon et alJ (|2000() . An order blocking filter (6601/15) 
isolated Ha emission at order 379 over a 9 arcminutc 
field of view. At this order, the spectral resolution is 
quite high (FWHM ~ 0.5A = 22.9 km s -1 ), but the 
wavelength satisfying the etalon interference condition 
varies radially across a given image (each of which 
corresponds to a given etalon spacing, h). Over the 
course of the observing run, the etalon spacing was 
changed to sample the full free spectral range (FSR = 
17. 4A) at each CCD pixel. It should be noted, however, 
that telescope pointing variations resulted in the FSR 
being incompletely sampled at some spatial locations, 
resulting in blank pixels later in the reduction process. 
The exposure time at each of the 72 etalon spacings was 
12 minutes. To avoid confusion from ghost images, the 
galaxy was placed away from the optical center (a faint 
ghos t image can be seen in the upper-left corner of Fig. 
lAlsK ). 

Bias frames, taken during the observing run, were av- 
eraged and subtracted from each image. Flat fields were 
averaged and applied to the images. A "white light" cube 
was also obtained. This reduction step is used to remove 
any wavelength-dependent flat field structure. However, 
the white light cube was found everywhere to vary by 
2 percent or less. A test application of the white light 
cube yielded little difference; therefore, this correction 
was deemed unnecessary. The individual images were 
arranged in order of etalon spacing and stacked into a 
data cube. To increase signal-to-noise, the images were 
spatially binned 2x2 in software. Cosmic ray removal 
was performed by hand using the IRAF 1 task CREDIT. 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
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The resulting cube contains planes of constant etalon 
spacing rather than wavelength. The surfaces of constant 
wavelength are paraboloids, centered on the optical axis 
and with a curvature constant K\ which is depende nt on 
instrumental par ameters (see iBland fc Tullvlll989l and 
our equation IA2|I . In addition, the position of the opti- 
cal axis on the detector is not known a priori; indeed, the 
position was found to vary over the course of the observ- 
ing run, distorting the shape of the constant-wavelength 
paraboloids. Moreover, telescope pointing variations re- 
sult in shifts of object locations with respect to both the 
detector and optical centers. The process of aligning the 
optical axis and transforming the surfaces of constant 
wavelength to planes is called the phase correction. Our 
procedure for subtracting night-sky line emission and ob- 
taining the wavelength solution (representing the bulk of 
the phase correction) is outlined in the Appendix. 

Relative intensity variations between the sky- 
subtracted images were measured and corrected for 
using field stars. The IRAF task GAUSS was used to 
convolve each image to a common beam corresponding 
to the worst seeing conditions from the observing run 
(3.55"). Because we are primarily looking for diffuse 
emission, the poor seeing is not harmful; rather, the 
faint diffuse structure is brought out by the convolution. 

Astrometric solutions for each image were calculated 
with the IRAF task CCMAP, using six relatively bright 
stars included in the HST Guide Star Catalog 2 . With 
these astrometric solutions, and the wavelength solution 
described in the Appendix, the phase correction was com- 
pleted by sorting each intensity value into the correct 
pixel in the final data cube. In cases where multiple in- 
tensity values were assigned to a single pixel, the average 
value was used. The pixels in the final data cube are 2" 
square, and the channel width is 11.428 km s -1 . 

With a FSR of only 17.4A, corresponding to ~800 km 
s _1 at Ha, the r otational velocity of the galaxy (198 km 
s -1 : Ilrwmlll994[) together with the broad emission pro- 
files (o-gas ~ 32.5 — 42.5 km s" 1 ; see EI3.2f) causes real 
Ha emission to appear in nearly all planes of the data 
cube. Thus, there are few continuum channels in our 
data set and the continuum must be removed locally. To 
remove the continuum emission, a median-filtering algo- 
rithm was implemented, using intensity values along the 
spectral axis at a particular spatial location and those 
along the four nearest-neighbor spectra. Taking the in- 
tensity values from all five spectra together, the median 
and standard deviation were calculated. Values differ- 
ing from the median value by more than two standard 
deviations were excluded (thus eliminating the part of 
the spectrum containing Ha emission), and new statis- 
tics were calculated. This procedure was repeated until 
no more statistical outliers were found, or for a maximum 
of five iterations. The resulting median was taken as a 
measure of the continuum level at that location, and was 
subtracted from the spectrum. 

To enhance faint emission far from the galaxy mid- 
plane, the data cube was further smoothed to an 8" 

for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

2 The Guide Star Catalog was produced at the Space Telescope 
Science Institute under U.S. Government grant. These data are 
based on photographic data obtained using the Oschin Schmidt 
Telescope on Palomar Mountain and the UK Schmidt Telescope. 



beam. Blank (unsampled) pixels in the final cube, caused 
by variations in telescope pointing, were replaced before 
smoothing by interpolating over adjacent pixels. In the 
region of the data cube containing emission from NGC 
5775, approximately 3 per cent of the pixels were initially 
blank. 

To obtain a rough intensity calibration, a major axis 
cut through the moment-0 map shown in Fig was 
compared to a majo r axis cut through the Ha image from 
Colli ns~et al.l l)2000(l . smoothed to the same resolution as 
in the moment-0 map. The noise in the channel maps was 
thus measured to be 1.38 x 10 -19 ergs cm -2 s _1 arcsec -2 
channel -1 . Assuming a gas temperature T=10 4 K, this 
corresponds to an emission measure (EM) per channel of 
0.0688 pc cm -6 channel -1 . 

Because the wavelength solution was based on night- 
sky emission lines rather than observations of a stan- 
dard wavelength calibration lamp (see the Appendix), 
the velocity scale was ve rified by com paring with exist- 
ing data. HI data fsee Ilrwinl [l994) and CO 2-1 data 
("see iLee et "all 12001)) were kindly provided by J. Irwin 
and S.-W. Lee for this purpose. The beam size of the 
HI data is 13.6" x 13.4" at position angle -33.7°, and 
the channel width is 41.67 km s _1 . The beam size of 
the CO 2-1 data is 21", and the channel width is 8.08 
km s . Overlays of major axis position-velocity (PV) 
diagrams were generated to compare both data sets to 
the Ha data. The HI PV diagram differs substantially 
from the Ha PV diagram in that it shows a slower rise 
in velocity at low R (see Fig. ^,). This is due to the HI 
being less centrally concentrated than the ionized g as , as 
found through modeling of the H I data cube by Irwin 
(1994). The HI diagram is therefore not optimal as a 
check on the velocity scale. On the other hand, the CO 
emission is expect ed to follow the ionized gas distribution 
more closely (e.g.. iRownd fc Youndll999HWong fc Blitz! 
2002), except where extinction may affect the Ha pro- 
files. Indeed, the shapes of the PV diagrams are more 
similar (see Fig. H 3 ). Based on the comparison between 
the CO and Ha PV diagrams (after converting the Ha 
velocities to the LSR frame), a constant 9 km s -1 off- 
set was added to the velocity axis of the Ha data cube. 
This correction is smaller than the channel width of the 
data cube. In all presentations of kinematic data in this 
paper, velocities are relative to the systemic velocity of 
NGC 5775 (cf. Table UJ. 

Moment maps were generated with the Groningcn Im- 
age Processing System (GIPSY) task MOMENTS by re- 
quiring that emission appears above the 3cr level in at 
least three velocity channels. The moment-0 (total in- 
tensity) map is displayed in Fig. [5^,, and the moment- 1 
(mean velocity) map is displayed in Fig. ^jp. 

3. ANALYSIS AND MODELING 

The primary goal of this work is to study the kine- 
matics of the ionized halo of NGC 5775. In this Section, 
we first present the velocity field and point out some in- 
teresting features and trends. Next, major axis rotation 
curves are obtained for both the ionized and molecular 
gas components from PV diagrams. Our analysis of the 
kinematic structure is then extended by modeling the 
halo component to search for a change in the rotation 
curve with height above the midplane. 
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Fig. 1. — (a) Overlay of HI (white contours) and Ha (black 
contours) major axis PV diagrams. HI contours run from 3.36 to 
43.7 K in increments of 6.72 K. Ha contours run from 8.28 X 10~ 18 
to 7.04 X 10 -17 erg cm -2 s _1 arcsec -2 channel -1 , in increments 
of 6.90 X 10 — 18 erg cm -2 s — 1 arcsec -2 channel -1 , (b) Overlay of 
CO 2—1 (white contours) and Ha (black contours) major axis PV 
diagrams. CO contours run from 0.03 to 0.33 K in increments of 
0.05 K. Ha contours are the same as in (a). Positive values of the 
major axis distance R correspond to the southeast side of the disk. 
The 9 km s — 1 offset described in the text has already been applied 
to the Ha velocity scale. 
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Fig. 2. — (a) Moment-0 map of the Fabry-Perot data cube. 
Contours illustrate the appearance of the brighter structures in the 
disk, (b) Moment-1 map of the Fabry-Perot data cube. Velocities 
are relative to the systemic velocity of 1681.1 km s — 1 . Contours 
run from -210 to 210 km s in increments of 30 km s — 1 (the 
SE and SW quadrants constitute the receding side). The systemic 
velocity contour is darkened. The morphological major and minor 
axes, determined using the center and position angle listed in Table 
HI are indicated by dashed lines. 



3.1. The velocity field 

Some insight into the global kinematic characteristics 
of the ionized component of NGC 5775 can be obtained 
by examining the velocity field. Figure |2b shows the 
moment-1 map calculated from the data cube. The fol- 
lowing features are observed: 

• In all four galaxy quadrants, mean velocities are 
seen to decrease with increasing distance from the 
major axis (up to ~ 10"). This initial decrease is 
attributed to viewing a rotating disk in projection, 
and is not indicative of a halo lag. To determine the 
location of the edge of the projected disk, we as- 
sume that the disk is axially symmetric, with a ra- 
dius equal to the semi- major axis length of « 100". 
Beyond this distance, the emission is observed to 
fall sharply, and we therefore assume that the disk 
has a sharp radial cutoff. 

• In the southwest and northeast quadrants (labeled 
SW and NE in Fig. for most contours the 
mean velocities continue to decrease beyond the 
extent of the projected disk. The EDIG in these 
regions is dominat ed by filamentary structures (cf. 
ICollins et a,l]l200?ih . 

• In the southeast quadrant (labeled SE in Fig. I^Jd), 



the velocity contours show an increase in mean ve- 
locity above the edge of the projected disk. The 
contours then appear to run roughly parallel to 
the minor axis. In the northwest quadrant (labeled 
NW in Fig. Eb), the behavior of the mean veloci- 
ties is confused by the presence of a closed contour 
structure, but it appears that some of the contours 
show a continued drop in mean velocity above the 
edge of the projected disk, while other contours run 
parallel to the minor axis. 

• The asymmetry of the velocity contours along the 
major axis, particularly in the southern (receding) 
side, is likely caused by significant dust extinction 
in the disk. The dust lane, which runs parallel to, 
but offset to the northeast from the major axis, 
preferentially obscures emission from gas more dis- 
tant from the observer along those lines of sight. 
Because the velocity projections along the line of 
sight are minimized at the near edge of the disk, 
we observe lower mean velocities to the northeast 
of the major axis than to the southwest, where ex- 
tinction is not as extreme. In other words, extinc- 
tion in the dust lane results in little or no emission 
reaching us from near the line of nodes. 

• The kinematic and morphological minor axes are 
no t paral l el in t he projected disk [also noted in H I 
hv lTrwil iTmol ]. We assume that the emission in 
the projected disk locations is in fact dominated 
by disk emission. The sense of the offset, together 
with the inclination angle of the galaxy and the as- 
sumption that the near side is tilted to the north- 
east with respect to the morphological major axis 
(based on the appearance of the dust lane) , is con- 
sistent with a radial inflow of gas in the disk, or 
may indicate the presence of a bar. If the former, 
there is no clear indication that it extends beyond 
the projected disk into the halo. NGC 5775 is clas- 
sified as barred in the Third Reference Catalogue o f 
Bright Galaxies (RC3; lde Vaucouleurs et alTl 991). 
although an examinat ion of the Two Micron All 
Sky Survey C2MASS: iKTemmann et a.niTflfll K- 
band survey image does not suggest bar signatures 
such as a box- or peanut-shaped bulge. 

Though the moment-1 map is useful for viewing overall 
trends, it cannot be used to derive rotation speeds for a 
nearly edge-on galaxy. Intensity-weighted velocities are 
affected by the rotation curve, varying projections of the 
rotation-velocity vector, and the gas distribution along 
any line of sight and thus are not indicative of rotation 
speed. A full consideration of the line profile shapes and 
three-dimensional density distribution (as well as an as- 
sessment of extinction) is necessary to derive robust ro- 
tation curves in the case of high inclination. 

3.2. Disk rotation curves 

Here, we wish to construct disk rotation curves from 
major axis PV diagrams generated from both the Ha and 
CO 2-1 data sets. Several methods have been suggested 
for recovering the rotation curve of a galaxy. In par- 
ticular, two useful methods in the stud y of edge-on spi- 
rals are the "envelope traci ng" method ijSancisi fc Alien! 
119791 iSofue fc Rubhll2001() and the "iteration method" 
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ijTakamiva fe Sofuell2002ft . The envelope tracing method 
calculates the rotation curve using the high- velocity edge 
of a PV diagram. The iteration method automates the 
procedure of generating a model galaxy (with speci- 
fied radial density profile and rotation curve) that best 
matches the major axis PV diagram of the observed 
galaxy. The benefit of using the iteration method is a 
more accurate recovery of the rotation curve at small 
radii, which, because of beam smearing and rapidly 
changing densities and velocities, is typically underes- 
timated by the envelope tracing method. However, our 
implementation of the iteration method employs the en- 
velope tracing method as part of its fitting procedure; 
hence, we describe here the specifics of the algorithms 
used in both methods, which we have written as MAT- 
LAB scripts. For reasons discussed in ij3.2.2l the iter- 
ation method was unable to accurately determine the 
major axis rotation curve for NGC 5775, but the results 
of the method were used as a starting point for a visual 
determination. 

Throughout this Section, the disk is assumed to con- 
sist of a series of concentric, axisymmetric rings, each 
of which is allowed to have a different gas density, rota- 
tional velocity, and velocity dispersion. We do not allow 
the disk to be warped, under the assumption that the 
Ha emission traces the star forming disk and not the 
outer parts. An examination of the channel maps did 
not reveal the signature of a warp along the line of sight 
(see, e.g.. iSwaters. Sancisi. fc van der Hulstl 11997ft. and 
the moment-0 map does not show evidence for a warp 
across the line of sight. Here, we briefly describe the 
specifics of the algorithms used. 

3.2.1. The envelope tracing method 

After the description of lSofue fc Rubinl 1)2001^ . we use 
the following algorithm. A PV diagram is generated, and 
the velocity dispersion (cr gas ) is estimated by measuring 
the width of the line profile at the highest radius. It 
is assumed that at the highest radius, only one ring of 
emission is being sampled, so velocity projections from 
additional rings do not broaden the profile. At each ra- 
dius, the value of the rotation curve is taken to be the 
velocity at the edge of the line profile, with a correction 
for the velocity dispersion, the channel width, and the 
inclination angle of the galaxy. 

Mathematically, the inte nsity of the envelope on the 
line profile is defined to be ijSofue fc RubinlfeOOlj) 

4nv = yVmax^+Ai (1) 

where I max is the maximum intensity in the line profile, 
I\ c is the lowest contour value, typically taken to be 3 
times the rms noise in the PV diagram, and 77 is a con- 
stant, normally taken to be in the range 0.2-0.5. On the 
side of the line profile farthest from the systemic veloc- 
ity, th e rotational velocity is taken to be (|Sofue fc Rubinl 

mm 



V r = (Ucnv - V sys ) / shl(i) - yj af nst + CT^, (2) 

where v cnv is the velocity at which / = I env , v sys is the 
systemic velocity, i is the inclination angle, and Uinst is 
the instrumental velocity resolution. 



3.2.2. The iteration method 

Our implementation of the iteration method makes use 
of the envelope tracing method and the suite of GIPSY 
tasks. Most importantly, the task GALMOD is used to 
generate a model data cube using user-specified radial 
and vertical density profiles, rotation curve, velocity dis- 
persion, and viewing angle. To create the initial model, 
we specify: the velocity dispersion, which is estimated 
in the same way as for the iteration method; the radial 
density profile, which is estimated using the GIPSY task 
RADIAL; and the initial guess for the rotation curve, 
which is estimated by using the envelope tracing method 
on the observed data set. The radial profile must be well 
matched to the data because both the density distribu- 
tion and the rotation curve affect the shape of the line 
profiles. The scale of the input radial density profile is 
fixed such that the signal-to-noise in the model is ap- 
proximately that measured in the data. To measure the 
noise in the model, the same inputs are used with two 
different random number seeds, and the standard devia- 
tion of the difference between the two runs of the model 
is calculated. 

Once the initial model has been ge nerated, we follow a 
proced ure based on that detailed bv iTakamiva fc Sofud 
( 2002). The envelope tracing method is used to derive 
a rotation curve from the initial model. That rotation 
curve is compared to the one derived from the data, and 
the differences between the two are added to the input 
rotation curve to generate a new model. This proce- 
dure is repeated until the difference between the rota- 
tion curves derived from the data and the most recent 
model satisfies a convergence criterion or until a max- 
imum number of iterations have been performed. The 
degree of success of this method is dependent on how 
well the other galaxy parameters (radial density profile, 
velocity dispersion, signal-to-noise ratio) are reproduced 
in the model. If these parameters are poorly specified, 
the iteration method will not be able to correctly repro- 
duce the observed PV diagram. 

Because RADIAL assumes a circularly symmetric den- 
sity distribution, it is unable to match the observed 
major-axis intensity profile, which shows clear signs of 
dumpiness (see, for example, the bright knot of emission 
at R — +90" in Fig. The deviations from axisym- 
metry are severe enough that the model PV diagrams 
are too poorly matched to the data to allow the iteration 
method to converge. We therefore abandon this method 
as a means of directly obtaining the rotation curve, and 
attempt to mitigate some of the effects of the clumpy 
density distribution. Better results are obtained (though 
less efficiently) by starting with the converged radial den- 
sity profile obtained with RADIAL and a rough rotation 
curve obtained with the iteration method after a few iter- 
ations. These results are used to generate a modeled ma- 
jor axis PV diagram. The observed and modeled PV di- 
agrams are compared by eye, and first the radial density 
profile is adjusted so that observed and modeled major 
axis intensity distributions are reasonably well matched 
except in regions of obvious dumpiness. Then, the rota- 
tion curve is modified until a good match is achieved. 

A major-axis rotation curve and a radial density profile 
have been obtained in this way from the Fabry-Perot 
data. Fig. [3] shows an overlay of the major axis PV 
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Fig. 3. — Overlay of Ha (white contours) and best-fit model 
(black contours) major axis PV diagrams. Contour levels for both 
are 10<t to lOlOcr in increments of 100o\ Positive values of the 
major axis distance R correspond to the southeast side of the disk. 




+100 -100 

R (arcsec) 

Fig. 5. — Overlay of CO 2-1 (white contours) and best-fit model 
(black contours) major axis PV diagrams. Contour levels for both 
are 3ct to 18a in increments of 3<r. Positive values of the major 
axis distance Ft correspond to the southeast side of the disk. 
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Fig. 4. — (a) Profile of density versus galactocentric radius R for 
the best-fit models - Ha (squares) and CO 2—1 (diamonds), (b) 
Rotation curve for the best-fit mo dels — Ha (sq uares) and CO 2-1 
(diamonds). The rotation curve of Irwin 11994) is plotted (dashed 
line) for reference. 



diagrams obtained from the data and the best-fit model. 
Fig. Qlshows the radial density profile and rotation curve 
in the best-fit model. To achieve the best match, the 
velocity dispersion in the model was set to 32.5 km s _1 
on the approaching side, and 42.5 km s _1 on the receding 
side. The kinematic center and systemic velocity listed 
in Tabled were used, and provided good agreement. 

To check that the radial density profile and rotation 
curve are reasonable and that our results are not biased 
by, for example, extinction, we have repeated this pro- 
cedure for the CO 2-1 data. Fig. shows an overlay of 
the major axis PV diagrams obtained from the data and 
the best-fit model. Fig. 0] shows the radial density pro- 
file and rotation curve in the best-fit model. A velocity 
dispersion of 15 km s _1 was required to obtain the best 
match. Good agreement was obtained without modify- 
ing either the systemic velocity or the kinematic center 
listed in Table In order to match the low signal-to- 
noise ratio of the CO data, the modeled emission is quite 
faint; the apparent asymmetry between the approaching 
and receding sides of the model contours in Fig. [21 is due 
to the noise in the model. 

Despite the fact that our Ha and CO radial profiles 
differ, the rotation curves are very similar. The mean 
difference between the Ha and CO rotation curves is 1.5 
per cent, the maximum difference is 29.5 per cent (at low 
i?, where the CO rotation curve rises more slowly than 



does the Ha rotation curve), and the rms difference is 
11.5 km s . This correspondence implies that our pro- 
cedure works well. We note that the match between the 
Ha PV diagram and the corresponding model is worse in 
the inner parts (R < 20") than for the CO PV diagram. 
In the former, the data appear to rise more slowly with 
a higher velocity dispersion. 

3.3. Halo rotation 

Having recovered a radial density profile and rotation 
curve for the major axis, we now move on to model- 
ing the halo Ha emission. PV diagrams are constructed 
along cuts parallel to the major axis, at various heights 
above the midplane z (in this paper, z is positive to the 
southwest of the major axis). To calculate model PV di- 
agrams, we have made a modification to the GIPSY task 
GALMOD to allow for a vertical gradient in azimuthal 
velocity: 



v(R,z) = 



v{R,z = Q)-f z \ 
v(R,z = 0) 



z\ — z ] for z > z 
for z < zq. 



(3) 



where dv /dz is a constant parameter in the model, with 
units [km s _1 arcsec -1 ]. Note that this model fixes the 
shape of the rotation curve as a function of height, and 
only changes the amplitude. We include a parameter 
(zq) that specifies the height at which t he rotational lag, 
begins in order to be consistent with iFraternali et alJ 
(2005), who find that the neutral halo of NGC 891 co- 
rotates with the disk up to z — 1.3 kpc, and shows a 
vertical gradient in azimuthal velocity above that height 
(although the authors state that it is possible that the 
observed corotation at this height is an effect of beam 
smearing). Limited information is available to constrain 
whether zo differs from zero in our model. For NGC 
5775, a height of 1.3 kpc corresponds to z = 11". How- 
ever, evidence discussed in this Section indicates that a 
gradient is already present at z = 10". We therefore 
set zq to be the (exponential) scale height of the galaxy 
model (of order 5"; see below), and show later that a 
choice of z = 0" does not change the derived gradient 
significantly. 



Kinematics of DIG in NGC 5775 



7 



3 

>> 

U 

(0 
*-> 

• i-H 

u 

cO 



o 



o 
o 
o 



: 1 1 1 1 1 


1 1 1 1 : 


: NE 


sw : 


■/ 
■ / 

■ / 

■ / 
■ / 
■ / 


\ ■ 
\ ■ 
\ ■ 


■ / 

_ ■ / 

1,1,1 


i , i 



-40 -20 20 
z (arcsec) 



40 



Fig. 6. — Comparison of the observed (squares) and modeled 
(solid line) vertical intensity profiles. Both profiles were obtained 
by averaging the vertical distribution of emission along the length 
of the disk. The modeled profile is normalized to the observed 
profile at z = 0. The exponential scale heights (as described in 
the text) are 5.8" on the southwest side and 4.5" on the northeast 
side. At z = —30", the model is significantly lower than the data; 
hence the 5.8" scale height is used to generate PV diagrams at that 
particular height. 



Such a vertical gradient in azimuthal velocity is consid- 
ered in H3.3.1I We also consider the effects of modifying 
the shape of the halo radial density profile f H3.3.2|l . the 
shape of the halo rotation curve f t!3.3.3|l . and the position 
angle and inclination of the halo f <l3.3.4[l . 

Assuming a circularly symmetric disk with a sharp op- 
tical cutoff radius of 100" (cf. Figure EH) and an in- 
clination angle of 86°, the edge of the projected disk 
along the m i nor ax is is located at a distance of 7". 
Collins"et al.l II2 000T) report that, based on the model- 
ing of B vun. Freeman, fc Kvlafis (1994), extinction ef- 
fects should be negligible at z > 600 pc = 5". Therefore, 
at minor axis distances greater than 10", we assume the 
effects of the projected disk and extinction are negligible. 
After smoothing the data to an 8" beam, reasonable PV 
diagrams can be constructed up to a height \z\ w 30", 
but beyond that point there is not enough reliably de- 
tected emission. Thus, the range of modeled heights is 
10" < \z\ < 30" (1.2 kpc < \z\ < 3.6 k pc). 

First, the model described in i|3.2.2l is considered with 
no vertical gradient in azimuthal velocity. We call this 
the cylindrical rotation (CR) model. The scale height is 
adjusted to match the minor axis intensity distributions, 
but the radial density distribution and rotation curve de- 
rived for the major axis are unchanged. The necessary 
scale heights are 5.8" for the southwest side of the disk 
(z > 0) and 4.5" for the northeast side of the disk (except 
for the z = —30" PV diagram, for which the 5.8" scale 
height was retained) . The observed vertical emission pro- 
file is compared to the model in Figure[f)| Fig. [7|displays 
overlays of the observed PV diagrams and the CR model 
PV diagrams for slices at heights z = ±10", ±20", ±30". 

Upon examining these diagrams, some clear features 
are: 



• A rotation velocity gradient appears to be neces- 
sary primarily on the northwest side (R < 0). We 
attempt to model this gradient in c!3.3.1l In some 
panels, a shift in systemic velocity seems to be nec- 
essary. Such an offset is also included in t !3.3.1l 
The kinematics appear different on the positive- z 
and negative- z sides of the halo. 

• The "knee" in the model PV diagrams, caused by 
a combination of high gas density and rotational 
velocity, appears to be unnecessary in some panels, 
but in most, it appears that this knee moves fur- 
ther from R = 0" with increasing z. The shift in 
this knee may be explained either by a change in 
the radial density profile (demonstrated in S13.3.2fl . 
perhaps due to a radial migration, or by a change in 
the shape of the rotation curve with height (demon- 
strated in H3.3.3JI . as might be expected as the in- 
fluence of the bulge potential diminishes. The knee 
could also be a signature of non-circular disk mo- 
tions associated with a bar potential, although as 
mentioned in 13.11 there is little evidence for a bar 
in 2MASS images. 

• In panel /, the bright knot of emission at R « 
±100" and v ~ 150 km s -1 correspo nds to a known 
llCollins et alJl200(l ILee et al.ll2001^ Ha extension 
of possibly tidal origin. This emission is clearly vis- 
ible in the moment-0 map in Fig. extending to 
the east away from the southeast edge of the disk. 
That the radial velocity of this complex appears to 
be significantly lower than that of the surrounding 
gas either means that the rotational velocity has 
dropped significantly or that the complex does not 
lie along the line of nodes. 

3.3.1. Azimuthal velocity gradient 

The slit spe c tra p resented by iRandl ((2000) and 
iTullmann et alJ l)2000[) were examined to provide a first 
estimate of the magnitude of a possible vertical gradi- 
ent in azimuthal velocity. Three slits are available, each 
measuring gas velocities on both sides of the midplane, 
so that the gradient may be estimated in six regions. 
Based solely on the mean velocities, the vertical gradient 
in azimuthal velocity was estimated to be ~ 1 — 2 km 
s -1 arcsec -1 . 

To test for the existence of a vertical gradient in az- 
imuthal velocity of the form described by equation 
we have added a 1 km s -1 arcsec -1 (or, equivalently, 
about 8 km s -1 kpc -1 ) gradient to the CR model. Fig- 
ure |H1 demonstrates the change in shape of the PV dia- 
grams when the gradient is included in the model. The 
flat part of the approaching side appears to be better 
matched. However, there is still an apparent problem 
matching the shape of the PV diagrams for R < 60" in 
the halo (see H3.3.2l and c!3.3.3l for possible explanations). 
It is not clear that adding the gradient to the receding 
side has improved the model. 

An adjustment of the systemic velocity appeared to be 
necessary as well in most panels. The CR model was 
further modified by adding a +10 km s -1 offset to the 
systemic velocity. This offset is already included in Fig. 
03 Though the offset is smaller than the channel width 
of the Fabry-Perot data cube, it seems to improve the 
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Fig. 7. — PV diagram overlays of H« (white contours) and CR model (black contours). PV diagrams are displayed at (a) z = +10", 
with contour levels 10<r to 210<r in increments of 40<r; (b) z = +20", with contours levels 10<r to 50ct in increments of 10cr; (c) z = +30", 
with contour levels 3<r to 18cr in increments of 3o"; (d) z = —10", with contour levels 10<t to 210<r in increments of 40<r; (e) z = —20", 
with contour levels 6cr to 36<r in increments of 6cr; (f) z = —30", with contour levels 3cr to 18ct in increments of 3<r. Positive values of 
galactocentric radius R and z correspond to the southeast side of the disk and the southwest side of the halo, respectively. 



match between model and data (for all regions of the 
halo). Note that this adjustment was not indicated by 
the major axis PV diagrams; it is only added to im- 
prove the agreement between data and model PV dia- 
grams in the halo. We also note that applying a shift in 
the systemic velocity without retaining the gradient in 
azimuthal velocity is not sufficient to match the shape of 
the PV diagrams in the halo. 

We have also generated models with different values 
of the azimuthal velocity gradient. A gradient as low as 
0.5 km s" 1 arcsec -1 is found to be too small to match 
the observed PV diagrams. Gradients with higher am- 
plitude (in particular, dv/dz — 2 km s" 1 arcsec" 1 ) are 
able to better match some individual regions (primarily 
the flat part of the PV diagrams on the approaching side 
for z — ±20"), but provide poor agreement at z — ±30". 
We conclude that dv/dz — 1 km s" 1 arcsec" 1 provides 
the best overall agreement with the data, but that higher 
values of the azimuthal velocity gradient may be appro- 
priate in localized regions of the halo. We note that for a 
model with this gradient but with zq = 0", the azimuthal 
velocities would be approximately 5 km s" 1 lower at all 
locations above one scale height, corresponding to a 3 
per cent decrease at z = 30" for a flat rotation curve of 
v(R,z= 0) = 200 km s" 1 . 



In Figs. and [SJ the "knee" in the modeled PV di- 
agrams is too prominent when compared with the data 
(except at z = +10", for R < 0"). We next make an 
attempt to match the shape of the PV diagrams in the 
halo, retaining the features of our best halo model thus 
far. The shape of the rotation curve is unchanged; only 
the radial density profile is modified (in the same way 
for each height z), as shown in Figure^, by the intro- 
duction of a central depression. We have not attempted 
to recover the actual radial density profile in the halo. 
Rather, we mean to illustrate that a radial density pro- 
file different from that obtained for the major axis bet- 
ter reproduces the changing shape of the PV diagrams 
with height. The PV diagram overlays are shown in Fig. 
1101 We note that the z — —30" panel suggests a cen- 
tral depression even more severe than we have modeled. 
Further evidence for a change in the radial gas distribu- 
tion with height will be discussed in cuts through the 
moment-0 map (shown in Figure IT?}) suggest that such a 
central hole may well be present in the halo, especially 
on the negative- z side. Moreover, comparisons with PV 
diagrams constructed from the ballistic base model (dis- 
cussed in see Figure ITS)) suggest that the changing 
shape of the observed PV diagrams is strongly affected 
by radial redistribution of matter in the halo. 



3.3.2. Modification of the halo radial density profile 



3.3.3. Modification of the halo rotation curve 
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Fig. 8. — PV diagram overlays of Ha {white contours) and model [black contours) where a vertical gradient in azimuthal velocity of 1 
km s — 1 arcsec -1 has been added to both sides of the disk, and the systemic velocity has been increased by 10 km s — . Contour levels are 
as in Fig. Q 
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Fig. 9. — (a) Comparison between the profile of density versus 
galactocentric radius R in the b ase model (dotted line) and the 
modification presented in (solid line), (b) Comparison be- 

tween the rotation curve i n the base model (dotted line) and the 
modification presented in (solid line). 



Because the radial density profile and rotation curve 
are partially coupled in the PV diagrams, the density 
profile may not be responsible for the change in shape 
of the halo PV diagrams. We next attempt to modify 



the form of the rotation curve such that the shape of the 
halo PV diagrams is better matched. In this case, the CR 
model radial density profile is unchanged. The 1 km s _1 
arcsec -1 vertical gradient in azimuthal velocity and +10 
km s _1 systemic velocity offset are also still included. 
We have not attempted to recover the exact shape of the 
rotation curve in the halo. Rather, we illustrate how a 
change in the shape of the rotation curve changes the 
shape of the PV diagrams in the halo. The modification 
to the rotation curve is shown in Figure^. The rotation 
speed rises roughly linearly for R < 50", as would be the 
case for a less centrally condensed potential. Figure ITT1 
shows the PV diagram overlays. The modification to the 
radial density profile f M.3.2|l appears to match the data 
somewhat better than changing the shape of the rotation 
curve, but we cannot exclude the latter possibility. 

3.3.4. Modification of halo position angle and 
inclination 

In Figure^, the appearance of the velocity contours in 
the SW and SE quadrants could be interpreted as result- 
ing from the halo being oriented at a slightly different po- 
sition angle relative to the disk. The effect is not clearly 
visible in the NE and NW q uadran t s, alth ough the mean 
velocities in both spectra of IRandl l|2000fl are consistent 
with this asymmetry in all four quadrants. Such a shift 
in the position angle in a lagging halo would make the 
contours run more perpendicular to the major axis in two 
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Fig. 10. — PV diagram overlays of Ho (white contours) and model (black contours) where the radial density profile has been varied as 
shown in Figure Contour levels are as in Fig. |7| 



opposite quadrants (e.g., NW and SE in Fig. ^Jd), and 
increase the angle between the contours and the minor 
axis in the other two quadrants (e.g., NE and SW in Fig. 
Eh)- Examination of the velocity contours indicates that 
the magnitude of such an offset in position angle neces- 
sary to produce the observed effect would be on the order 
of 20° , but it should be noted that dumpiness in the halo 
gas distribution and peculiar velocities may confuse the 
situation. 

On a similar note, if we allow the inclination of the 
halo to vary relative to that of the disk, it is possible 
that an apparent lag in halo rotation could be simply an 
indication that the halo is viewed from a more face-on 
perspective than the disk. If we assume a disk rotation 
speed of 200 km s -1 and a vertical gradient in azimuthal 
velocity of 1 km s -1 arcsec -1 beginning at a height of 5", 
then at a height z = 30", the azimuthal velocities in the 
halo would be about 175 km s _1 . To mimic this effect 
with a variation in inclination angle, an offset between 
the disk and halo of about 30° would be required. We 
consider such an offset very unlikely. An offset in the 
position angle or inclination of the halo relative to the 
orientation of the disk could occur if the extraplanar gas 
were accreted. This possibility is interesting given the 
interaction between NGC 5775 and its neighbor, but it 
is uncertain how long such a configuration might last. 



4. THE BALLISTIC MODEL 



We next utilize the ballistic model of CBR to attempt 
to understand what vertical gradient in azimuthal veloc- 
ity would be expected if the disk-halo flow in NGC 5775 
is purely ballistic in nature. The full details of the model 
are described by CBR, who compared mean velocities 
from this model to those obtained from slit spectra for 
NGC 5775 and NGC 891. Basically, clouds are launched 
from the disk with an initial velocity selected from a con- 
stant probability distribution between zero and a maxi- 
mum "kick velocity", along a unit vector normal to 
the midplane for models considered here. The clouds are 
then allowed to orbi t ballistically i n the galactic gravi- 
tational potential of iWolfire et alJ ((1995). Whenever a 
cloud leaves the simulation or returns to the disk, it is 
replaced by another cloud. Initial locations of the clouds 
are randomly selected from a distribution with a radial 
exponential scale length i?o and a vertical gaussian scale 
height zq — 0.2 kpc. Other model parameters described 
by CBR are not used here. The clouds do not interact 
with each other, and are assumed to have constant tem- 
perature, density, and size (and therefore equal Ha in- 
tensities) throughout the course of their orbits. After the 
simulation is run for 1 Gyr, at which time the system has 
reached a condition of steady state, positions and veloc- 
ities of each cloud are extracted from the model. These 
outputs can be examined directly or used to generate an 
artificial data cube at any inclination, from which, e.g. 
PV diagrams or runs of mean velocity versus projected 
height above the midplane can be created. 
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Fig. 11. — PV diagram overlays of Ha (white contours) and model (black contours) where the shape of the rotation curve has been 
varied as shown in Figure |5J). Contour levels are as in Fig. [7] 



TABLE 2 
Ballistic Base Model 
Characteristics for 
NGC 5775 



Parameter 


Value 




6 kpc 


20 


200 pc 


v k 


160 km s" 1 




198 km s" 1 



The most influential parameter in the model is the ra- 
tio of the maximum kick velocity to the rotational veloc- 
ity of the disk, Vk/V c . For a given value of the circular 
velocity and initial R, increasing kick velocities result in 
more radial movement, larger maximum height of the or- 
bit, and increased drop in azimuthal velocity at the peak 
height of the orbit. The circular velocity is determined 
observationally, and the kick velocity is set such that the 
resulting scale height of cloud density matches the scale 
height of EDIG emission (determined by CBR) . This crit- 
ical parameter is thus reasonably well constrained by ob- 
servations. The other parameters were found to have 
little effect on the model outputs. The characteristics 
of the so-called base model for NGC 5775 from CBR are 
summarized in Tabled With these parameters set in the 
model, CBR compared the mean velocities obtained from 



the model viewed at i = 86° with those measured along 
their two slits for NGC 5775. They found that the mean 
velocities from the model were roughly the same as the 
measured ones, but the model could not reproduce the 
observed mean velocities at the largest heights, which are 
seen to approach the systemic velocity (see their Figure 
8). 

We have performed an analysis of the variation in az- 
imuthal velocity as a function of z in the ballistic model. 
This analysis shows that in the model that best matches 
the mean velocities from the slit spectra of CBR, the ver- 
tical gradient in azimuthal velocity is shallower than the 
corresponding variation in mean velocity. The reason for 
this discrepancy is the large-scale radial redistribution of 
clouds in the ballistic model. Figure IT21 displays contour 
plots of cloud density in the ballistic model as a function 
of R and z, and Figure ^5] shows the azimuthal veloc- 
ity curves of the ballistic model clouds as a function of 
height. Most of the clouds at high z are also found at 
large R. Therefore, most of the decrease in mean veloc- 
ity in the model is caused by velocity projection, and the 
magnitude of the vertical azimuthal velocity gradient is 
of lesser importance in setting the mean velocity gradi- 
ent. For this reason, the results of the ballistic model can- 
not be relied upon to explain variations in mean velocity 
unless the observed radial density profiles are found to be 
similar to those in the model. In the present study, our 
spectral resolution is sufficient to allow the density distri- 
bution to be modeled, thus allowing azimuthal velocities 



12 



Heald et al. 



$ 6 



: (a) 
_ v z > 


■ 

Ja " ■ 






_ v z < 












(c) 
-All v s 


ff^'iTiy- 











R (kpc) 

Fig. 12. — Contour plots of cloud density as a function of galac- 
tocentric radius R and z in the ballistic "base model" (see text), 
for clouds (a) moving up, (b) moving down, and (c) all clouds. 
Contours levels correspond to 80, 400, 880, 1520, 2320, and 3280 
clouds kpc -2 in (a) and (b). In (c), contour levels correspond to 
160, 800, 1760, 3040, 4640, and 6560 clouds kpc" 2 . 
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Fig. 13. — Plots of average azimuthal velocity versus galactocen- 
tric radius R in the ballistic base model, for clouds (a) moving up, 
(b) moving down, and (c) all clouds. Points are plotted in (a) for 
2 = kpc (diamonds), z = 2 kpc (circles), z = 4 kpc (crosses), 
z = 6 kpc (squares), 2 = 8 kpc (triangles), and z = 10 kpc (plus 
signs). The same symbols are plotted in (b) and (c), except that 
the 2 = 8 kpc and z = 10 kpc points are left out for clarity (the 
azimuthal velocities at those heights are very similar to the ones 
at 2 = 6 kpc). 



to be measured as a function of z and compared directly 
with azimuthal velocities from the ballistic model. 

As an aside, we point out how the azimuthal velocities 
to be extracted from the model may depend on the phys- 
ical nature of the disk-halo flow. In FigureEl we display 
meridional plots of cloud orbits in the ballistic model at 
various initial radii. The ballistic model predicts that 
the radial motion during the majority of a cloud's or- 
bit is radially outward. Only at the highest initial radii 
and kick velocities do the ends of the orbits show a radi- 



ally inward motion. Because the radius of a given cloud 
is nearly always increasing, conservation of angular mo- 
mentum dictates that the azimuthal velocity of the same 
cloud is nearly always decreasing. The azimuthal veloci- 
ties of upward-moving clouds are thus always higher than 
those of downward-moving clouds. The vertical gradi- 
ent in azimuthal velocity extracted from the model will 
therefore be dependent on the assumption of the dynam- 
ics of the disk-halo flow. In a scenario where the gas be- 
gins in the ionized state, cools and condenses into neutral 
clouds somewhere near the top of the clouds' orbits, the 
azimuthal velocities will be relatively high. On the other 
hand, if the flow begins as hot gas and cools to a warm 
ionized gas for the downward portion of the flow, the az- 
imuthal velocities will be lower. In the latter case, the 
assumption of ballistic motion is likely violated for the 
upward portion of the flow, and simply associating the 
warm ionized gas with the downward moving clouds in 
our model may not be accurate. The fact that the struc- 
ture of the EDIG emission resembles shells and filaments 
implies that at least part of the ionized distribution is 
upward moving, but additional evidence is required to 
determine how much of the flow is ionized. 

Clearly, the vertical gradient in azimuthal velocity in 
Figure EH is n °t as simple as that described by equation 

13 which assumes a constant rotation curve shape as a 
function of height. Nevertheless, an approximate value 
of dv/dz can be obtained by inspecting the plots. In the 
interest of comparison with the velocity gradient mod- 
eled from the data, we only report here the gradient seen 
in the ballistic model for kpc < z < 4 kpc, which is a 
slightly larger vertical range than we have modeled in this 
work (Sj^J, and for 5 kpc < R < 12 kpc. (No clouds in the 
model reach z = 4 kpc for R < 5 kpc, so azimuthal veloc- 
ities cannot be extracted at those locations. The upper 
limit, 12 kpc, corresponds approximately to the optical 
radius of NGC 5775.) For the upward-moving clouds, 
the average gradient over the indicated ranges of R and 
z is approximately 2.4 km s -1 kpc -1 ; for the downward- 
moving clouds, 4.3 km s -1 kpc -1 ; for all clouds consid- 
ered together, 2.6 km s -1 kpc -1 . Even the highest value, 
4.3 km s -1 kpc -1 (appropriate if all of the observed Ha 
emission is from clouds returning to the disk) , is approx- 
imately a factor of 2 lower than the vertical gradient 
modeled from the data [about 8 km s -1 kpc -1 ; recall 
that a galaxy model using dv/dz ss 4 km s -1 kpc -1 was 
examined and rejected (cf. H3.3.1J) ]. 

The base model of CBR includes a ratio of Vk/V c — 
0.81 in order to match the observed scale height of EDIG 
emission. To understand how the maximum kick ve- 
locity affects the vertical gradient in azimuthal velocity, 
and thereby to examine whether a model with a differ- 
ent value of Vfe might better match the gradient esti- 
mated from the data, we have generated models with 

14 = 100 km s -1 {Vk/V c = 0.51) and V k = 220 km 
s -1 (Vk/V c = 1.11). In those cases, the gradients in 
azimuthal velocity were found to be approximately 1.9 
km s -1 kpc -1 and 6.8 km s -1 kpc -1 , respectively, for 
downward-moving clouds only. The approximately lin- 
ear relationship thus determined between the gradient in 
azimuthal velocity and Vk/V c implies that a maximum 
kick velocity of 255 km s -1 (Vk/V c — 1.29) is required to 
match the observed gradient even when only downward- 
moving clouds are considered. Unfortunately, this value 
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Fig. 14. — Meridional plots of clouds in the ballistic base model with kick velocities equal to the maximum value, = 160 km s 1 , 
starting at galactocentric radii (a) R = 4 kpc; (b) R = 8 kpc; (c) R = 12 kpc; and (d) R = 16 kpc. Cloud positions are plotted at 20 Myr 
intervals. Azimuthal velocities (in km s _1 ) are noted above the cloud's plotted position for most time steps (some omissions are made for 
clarity) . 



of Vk/V c yields a DIG layer with a vertical exponential 
scale height of approximately 6.3 kpc, or about a factor 
of 3 higher than the values derived from the two slit spec- 
tra used by CBR of 2.1 — 2.2 kpc. We conclude that it is 
unlikely that any value of 14 will result in a model that 
reproduces the observed gradient in azimuthal velocity. 

CBR also consider two galactic potential models (2 and 
2i) from lDehnen fc Binnevl 1)19981) . They find that model 
2i, which has an oblate, flattened halo (axial ratio q = 0.3 
rathe r than q = 0.8 as in model 2; see lDehnen fc Binnevl 
1998), produces a steeper change in radial migration with 
the ratio Vk/V c compared to their base model. We might 
therefore expect a steeper gradient in azimuthal veloci- 
ties with height if we use this potential. Referring to 
Figures 4 and 5 of CBR, we select Vk = 200 km s _1 
for model 2i, which should give roughly the same ra- 



dial redistribution as the base model with Vk = 255 km 
s , for radii R > 5 kpc. The gradients in azimuthal 
velocity generated by this model are approximately 3.2 
km s _1 kpc -1 for downward-moving clouds, 1.0 km s _1 
kpc^ 1 for upward-moving clouds, and 1.6 km s" 1 kpc -1 
for all clouds considered together. This model produces 
a disk with an average vertical scale height of about 3.5 
kpc (considerably higher than the observed scale height, 
2.1 — 2.2 kpc). Increasing the maximum kick velocity 
still higher to increase the azimuthal velocity gradient 
will only make the scale height larger. Thus, a flattened 
dark halo is not likely to be more successful in explaining 
the observed gradient in azimuthal velocity. 

Figure IT31 demonstrates how the gradient in mean ve- 
locity in the model is affected not only by the gradient in 
azimuthal velocity but also by the radial redistribution 



14 



Heald et al. 




(•>) 0- 

R - 12 kpc 



5 10 

z (kpc) 



Fig. 15. — Plots of mean velocity as a function of height above 
the midplane (z) in the ballistic base model for clouds at major axis 
distances (a) R = 4 kpc and (b) R = 12 kpc. Mean velocities are 
shown for clouds moving up (squares), moving down (diamonds) , 
and all clouds (solid lines). All mean velocities were calculated 
assuming an inclination angle of 86° . 
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of clouds. The gradient in mean velocity versus height 
in the model is shown for R — 4 kpc and R = 12 kpc. 
Fig. 115b (R, = 4 kpc) shows a much steeper gradient in 
mean velocity with height than Fig. 115b (R — 12 kpc). 
The average gradients over the range kpc < z < 4 kpc 
are approximately 14 km s -1 kpc -1 in the former case 
and 6 km s^ 1 kpc -1 in the latter, despite the gradient 
in azimuthal velocity being roughly the same at these 
two radii. In Fig. ^1 the cloud distribution reaches 
3 — 4 times higher at R = 12 kpc than at R = 4 kpc. 
Thus, the only clouds encountered along a line of sight 
at low radius and high z are concentrated at the edges 
of the disk, whereas a line of sight at higher radius and 
the same z encounters more clouds closer to the line of 
nodes, keeping the mean velocity closer to the actual ro- 
tation speed at that height. At both radii, the gradient 
in mean velocity is much larger than that of azimuthal 
velocity because of radial redistribution. 

PV diagrams constructed from the output of the bal- 
listic base model and displayed in Figure fTHl clearly illus- 
trate the points previously discussed and may be directly 
compared to the data. Along the major axis (Fig. fTBkl. 
the appearance of the "knee" in the ballistic model PV 
diagram is roughly matched to that in the data. We 
note that the radial density profile and major axis rota- 
tion curve in the ballistic model have not been adjusted 
to the extent described in an d that the ballistic model 
does not simulate noise, so signal-to-noise contours may 
not be plotted and directly compared to the data. As in 
the observations, the "knee" in the ballistic model PV 
diagrams moves radially outward with increasing height 
above the midplane. Simultaneously, the velocity profile 
peak at a given R shifts away from the local azimuthal ve- 
locity and closer to the systemic velocity. Both effects are 
caused primarily by the radial redistribution of clouds in 
the halo, and appear more striking in the ballistic model 
than in the data (see also Figure ITTjl . The PV diagrams 
in Fig. 1 161 make clear that, at least in the ballistic model, 
most of the changes are caused by radial motions reshap- 
ing the velocity profiles. A gradient in azimuthal velocity 
is present in the model, and, though clearly evident by 
inspecting the shape of the PV diagrams, produces rel- 
atively minor modifications. That the effects described 



Fig. 16. — PV diagrams constructed from the output of the bal- 
listic base model (black contours), which has been projected to the 
distance of NGC 5775, smoothed to an 8" beam, and viewed at 
i = 86°. PV diagrams were created at heights of (a) z = 0", (b) 
z = 10", (c) z = 20", and (d) z = 30". Contour levels correspond 
to 0.5, 2, 3.5, 5, 6.5, and 8 clouds arcsec -2 channel -1 in (a); 0.5, 
1.5, 2.5, 3.5, 4.5, and 5.5 clouds arcsec -2 channel -1 in (b); 0.5, 1, 
1.5, 2, and 2.5 clouds arcsec -2 channel -1 in (c); 0.5, 0.75, and 1 
clouds arcsec -2 channel -1 in (d). Also plotted in each frame are 
PV diagrams constructed from the Ha data cube (white contours) 
on the positive-z (southwest) side of the halo. Contour levels cor- 
respond to 20<7 to 520<t in increments of 50a in (a); 20<r to 200cr 
in increments of 30<r in (b); 5cr to 35<r in increments of 6<r in (c); 
2<r to llcr in increments of 3cr in (d). The channel width in the 
ballistic model is the same as that of the Ho data cube, 11.428 km 
s . Positive values of the major axis distance R correspond to 
the southeast side of the disk. 



above are also seen in the observed PV diagrams pro- 
vides further evidence for radial redistribution in NGC 
5775. 

Whether the large-scale radial redistribution of gas 
predicted by the ballistic model is actually observed in 
NGC 5775 is a very important question. Such a redis- 
tribution may be apparent in cuts taken parallel to the 
major axis through moment-0 maps at various heights 
above the plane. Figure IT7I shows such cuts for momcnt- 
maps made from the base ballistic model and the Ha 
data cube. At most heights, it appears that the data 
cuts may follow the shape of the model cuts fairly well, 
apart from the complication of filamentary structures. 
This correspondence is better for z < 0". At heights 
greater than 30", the moment-0 cuts from the ballistic 
model show a much clearer signature of radial redistri- 
bution, but our data are not sensitive enough to make 
a comparison at those heights. We conclude, then, that 
there is moderate morphological evidence for radial re- 
distribution at the level predicted by the ballistic model. 

CBR have previously compared the results of this bal- 
listic model to mean velocities obtained with slit spectra 
for NGC 891 and NGC 5775. In NGC 891, the authors 
find that the observed mean velocities drop more slowly 
as a function of height than the mean velocities derived 
from the ballistic model. The authors concluded from 
this that (magneto-) hydrodynamical effects are proba- 
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Fig. 17. — Comparison between moment-0 cuts generated from 
the Ha data (squares) and the ballistic base model (solid lines), at 
heights of (a) z = -30", (b) z = -20", (c) z = -10", (d) z = 0", 
(e) 2 = +10", (f) z = +20", and (g) z = +30". In each panel, the 
mean of the model profile is scaled to the mean of the data. The 
major axis distance R is positive on the receding (southeast) side 
of the disk, and z is positive to the southwest. 



bly at work. For instance, the radial migration could be 
modified by a gas pressure gradient in the halo. However 
if, as seems reasonable, halo gas pressure decreases with 
radius, the discrepancy would be worse because outward 
radial migration would be larger and modeled mean ve- 
locities even closer to systemic. It is possible that one 
major reason for the difference between mean velocities 
from the data and the model is that the radial density 
distribution in the halo of NGC 891 is more centrally 
concentrated than the prediction of the ballistic model. 
That outward radial migration o ccurs a t some level in the 
halo is indicated by Figure 11 o However, 
to properly examine halo rotation, azimuthal velocities 
must be estimated from data and the density distribu- 
tion modeled, as in the present study. In a forthcoming 
paper, we attempt such an analysis from high spectral 
resolution data of the DIG halo of NGC 891 with two- 
dimensional spatial coverage. 

In the case of NGC 5775, CBR find a better overall 
agreement between the mean velocities obtained from 
the slit data and those from the ballistic model. How- 
ever, the ballistic model does not reproduce all of the ob- 
served features. Of particular interest with respect to the 
present study, the modeled mean velocities in the range 
z = — 4 kpc (the same range considered here) show a 
markedly steeper gradient than the data. Mean veloci- 
ties calculated from the Fabry-Perot data cube show a 
similar lack of agreement with those from the ballistic 
model. In both cases, small-scale velocity and density 
variations are likely significantly affecting the observed 
mean velocities. 

5. CONCLUSIONS 

We have obtained Fabry-Perot spectra of the Ha emis- 
sion line in the nearly edge-on spiral galaxy NGC 5775. 
Major axis radial density profiles and rotation curves 



were obtained for the ionized and molecular components 
of this galaxy. The observations have also allowed us to 
examine the azimuthal velocity variation as a function of 
height in the halo. We have found that a vertical gra- 
dient in the azimuthal velocity with a magnitude 1 km 
s _1 arcsec -1 , or about 8 km s -1 kpc -1 , is able to re- 
produce the gross features of PV diagrams constructed 
parallel to the major axis (but a larger gradient may be 
appropriate in localized regions of the halo) . Such a gra- 
dient is primarily indicated for the approaching s ide of 
the ga laxy, though mean velocities from slit 2 of iRandl 
( 2000) suggest the presence of a gradient on the reced- 
ing side as well (at least at the highest z, but recall that 
the radial density profile was not taken into account in 
that analysis). The magnitude of the gradient should be 
considered approximate because of uncertainties in the 
radial density profile and rotation curve adopted (espe- 
cially for the halo gas) in this study. One should also be 
aware that the interaction with NGC 5774 may lead to 
significant deviations from axisymmetry in the halo, as 
has been assumed in the modeling here. Nevertheless, 
it is apparent that a non-zero gradient is present in this 
galaxy. Comparisons between PV diagrams constructed 
from the data and from galaxy models suggest either a 
radial redistribution of gas in the halo or a shallower rise 
in the rotation curve than is observed along the major 
axis. Further evidence for the former possibility is pro- 
vided by total intensity profiles constructed at various 
heights parallel to the major axis. 

The ballistic model of CBR has been analyzed in more 
detail. We have found that while azimuthal velocities de- 
crease gently with increasing z, the corresponding mean 
velocities decrease more steeply (particularly at lower ra- 
dius). This steep gradient in mean velocity is due to 
radial outflow of gas in the ballistic model. This re- 
sult emphasizes the importance of using caution when 
interpreting mean velocities in edge-on or nearly edge-on 
galaxies. If the radial density profile of radiating gas is 
not well understood, the use of mean velocities as indi- 
cators of rotation speed can be highly misleading. 

The decrease in azimuthal velocity with height in the 
ballistic model is, in fact, shallower than that which has 
been inferred from the data, suggesting that additional 
mechanisms are important. We note here some effects 
which may be at work in the halo of NGC 5775. In 
a picture of a fluid disk in hydrosta tic equilibrium, as 
described by, e.g., iBeniamiil 1)200(1) . the steeper gra- 
dient inferred from the data may imply pressure de- 
clining with radius in the halo, which is not unrea- 
sonable to expect. By considerin g the baroclinic solu- 
tions to stationary hydrodynamics. fBarnabe et alJ (2005) 
were successful in generating a model in agreement with 
the vertical gradient in the rotation curve of NGC 891 
l)Fraternali et alJ 120 05): it would be interesting to test 
whether this method is able to reproduce the gradient 
measured in this work as well. A completely different 
physical pictu re which must also be co nsidered is that of 
gas accretion. IKaufmann et"aTI |2005) were also able to 
reproduce the lag in the halo of NGC 891 with SPH sim- 
ulations of infalling multiphase gas. Although we have 
considered the halo in this actively star forming galaxy 
to be star formation driven, the fact that NGC 5775 is 
interacting with its companion may suggest that such 
considerations are relevant here. 
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APPENDIX 

SKY SUBTRACTION AND WAVELENGTH SOLUTION DETERMINATION FOR FABRY-PEROT 

OBSERVATIONS 

The process of converting a Fabry-Perot data cube containing planes of constant etalon spacing (h) to a cube with 
planes of constant wavelength is called the phase correction. To complete the phase correction, the sky coordinates 
and wavelength of each pixel in the original cube must be specified. In this Appendix, we describe our method for 
subtracting the night sky line emission from the data, and determining the wavelength solution (the bulk of the 
phase correction process) . Typically, arc-lamp comparison spectra are observed to facilitate obtaining the wavelength 
solution. Unfortunately, an arc-lamp cube was not obtained during our observing run. Therefore, night-sky emission 
lines were used for the wavelength calibration. Because the wavelength satisfying the etalon interference condition 
varies radially across each image, the night-sky emission lines appear as rings in the images (see Fig. IA18h), w ith radii 
dictated by the wavelength of the line, the etalon spacing, and the curvature constant K\ (cf. equation IA2JI . Thus, 
by measuring the properties of a ring (center and radius) in each frame so that it may be subtracted, the wavelength 
solution is also obtained. A first attempt at findi ng the optical center of each image was made by fitting a circle to 
three points on an individual ring, as described bv lBland fc Tullvl (^989). However, the rings appear at about the 3.5(7 
level in the images, and because the widths of the ring pro files increase with decreasing ring radius it w as often difficult 
to select appropriate points for a good fit. As noted bv lJones. Shopbell. fc Bland-Hawthornl l)2002j) . the azimuthal 
symmetry of the ring images may be exploited to find the optical center and subtract the sky rings. The location of 
the optical center in each plane was found with the following algorithm: 

1. First, the location of the optical center is approximated by fitting a circle to three points on an individual ring. 
This step yields a point close to the true optical center. 

2. Next, object signal is masked, and an azimuthal average is performed about the fiducial ce nter selected in step 
(1). This results in a radial profile with clear ring signatures. An example is shown in Fig. IA18b . 

3. Next, a non-linear least-squares algorithm is used to fit a gaussian to the most prominent ring profile. 

4. The fiducial center point is then varied over a search grid. At each grid point, a new azimuthal average is 
calculated, and a gaussian is fit to the new ring profile. At the end of this step, the width of the line profile has 
been calculated at each fiducial optical center location. 

5. Because of the azimuthal symmetry, the profile width is minimized at the true optical center. The grid of line 
widths obtained in step (4) is therefore interpolated, and the location of the minimum on the interpolated grid 
is taken to be the true optical center. 

The azimuthal profile corresponding to our choice of optical center is used to form a two-dimensional image of the 
sky rings (Fig. IA18h ). which is then subtracted from the original, resulting in a sky-subtracted image (Fig. IA18H ). 
Azimuthal variations in ring intensity and any non-circularity of the rings will result in incomplete subtraction of the 
sky emission. The residual emission was minimized by limiting the image area to a 4.6 arcminute square, inside of 
which the azimuthal variations were greatly reduced. The calculation of the azimuthal average was limited to this area 
in each step, and the sky-subtracted image was cropped to the same region. 

The ring fitting algorithm also yields the radii of the sky rings in each imag e. The two brightest s ky lines were 
identified as OH A6596.64 and the unresolved pair OH AA6603.99, 6604.28 (cf. lOsterbrock et al.lll996Tl . We used a 
central wavelength of 6604. 12A for the unresolved pair. For each image, the sky-line wavelengths A and radii r are 
inserted in 

A = A + K\r 2 (Al) 

to determine Ao, the wavelength at the optical center. The curvature constant, K\, is dependent solely on instrumental 
parameters: 

K X = (n/2)AA (p2//2 am ) ; (A2) 

where n is the etalo n order, AAn is the on -axis free spectral range, is the CCD pixel si ze, a nd / ca m is the focal length 
of the camera lens l|Bland fc Tullvl lT989). Once Ao is known for each image, equation IA1I can be used to determine 
the value of A at every pixel (using its distance, r, from the optical center). The value of K\ [actually, its counterpart 
expressed in terms of etalon spacing, Kh — (n/2)Aho(pf l /f 2 am )], was verified by inserting the nominal values for the 
instrumental parameters and comparing to the curvature value obtained by fitting a simple parabola to a plot of 
night-sky emission ring radius versus etalon spacing. The values of Kh obtained with these two methods were found 
to differ by only 0.2 per cent; therefore, the nominal instrumental parameters were used to specify K\. 



Kinematics of DIG in NGC 5775 



17 



400 



300 



200 




200 



300 



400 



400 



300 



200 






o 

o 



CO 

a 




100 200 300 

Image radius (pixels) 



300 



200 



100 




100 



200 



300 



200 



300 



400 



Fig. A18. — (a) An original image at an individual etalon spacing. Note the prominent sky rings, and the faint ghost image in the 
upper-left corner, (b) Azimuthal average of the image displayed in a, with the object pixels masked out. (c) The azimuthal average has 
been used to generate an image of the sky rings, (d) To generate this image, the image displayed in c has been subtracted from that in o. 
The same grayscale values have been used to display each image. 



Once the value of A is known for each pixel, and the astrometric solutions have been obtained for each plane of the 
data cube (cf. Sj2Jl, the phase correction may be completed. 



REFERENCES 



Barnabe, M., Ciotti, L., Frat ernali, F., & Sancisi, R. 2005, A&A, 

in press (astro-ph/0509021) 
Benjamin, R. A. 2000, Rev. Mexicana Astron. Astrofis. Ser. Conf. 

9, Astrophysical Plasmas: Codes, Models, and Observations, ed. 

J. Arthur, N. Brickhouse, & J. Franco (Mexico, DF: Inst. Asron., 

UN AM), 256 
Bland, J. & Tully, R. B. 1989, AJ, 98, 723 
Bregman, J. N. 1980, ApJ, 236, 577 
Bregman, J. N. & Pildis, R. A. 1994, ApJ, 420, 570 
Byun, Y. I., Freeman, K. C, & Kylafis, N. D. 1994, ApJ, 432, 114 
Collins, J. A., Rand, R. J., Duric, N., & Walterbos, R. A. M. 2000, 

ApJ, 536, 645 



Collins, J. A., Benjamin, R. A., & Rand, R. J. 2002, ApJ, 578, 98 
(CBR) 

de Vaucouleurs, C, de Vaucouleurs, A., Corwin, H. C, Jr., Buta, 
R. J., Paturel, C, &; Fouque, R. 1991, Third Reference Catalogue 
of Bright Galaxies (New York: Springer) 

Dehnen, W. & Binney, J. 1998, MNRAS, 294, 429 

Dettmar, R.-J. 1990, A&A, 232, L15 

Dettmar, R.-J. 2005, in ASP Conf. Proc. 331, Extra-Planar Gas, 

ed. R. Braun (San Francisco: ASP), 155 
Fraternali, F., Oosterloo, T., Sancisi, R., & van Moorsel, G. 2001, 

ApJ, 562, L47 

Fraternali, F., Cappi, M., Sancisi, R., k. Oosterloo, T. 2002, ApJ, 
578, 109 



18 



Heald et al. 



Fraternali, F., Oosterloo, T., & Sancisi, R. 2004, A&A, 424, 485 
Fraternali, F., Oosterloo, T. A., Sancisi, R., & Swaters, R. 2005, 

in ASP Conf. Proc. 331, Extra-Planar Gas, ed. R. Braun (San 

Francisco: ASP), 239 
Gordon, S., Koribalski, B., Houghton, S. & Jones, K. 2000, 

MNRAS, 315, 248 
Hoopes, C. G., Walterbos, R. A. M., & Rand, R. J. 1999, ApJ, 522, 

669 

Irwin, J. A. 1994, ApJ, 429, 618 

Jones, D. H., Shopbell, P. L. & Bland-Hawthorn, J. 2002, MNRAS, 
329, 759 

Kaufmann, T., Mayer, L., Wadsley, J., Stadel, J., & Moore, B. 

astro-ph/0507296 
Kleinmann, S. (j., Lysaght, M. G., Pughe, W. L., Schneider, S. E., 

Skrutskie, M. F., Weinberg, M. D., Price, S. D., Matthews, K., 

Soifer, B. T., & Huchra, J. P. 1994, Ap&SS, 217, 11 
Lee, S.-W., Irwin, J. A., Dettmar, R.-J., Cunningham, C. T., Golla, 

G., & Wang, Q. D. 2001, A&A, 377, 759 
Miller, S. T. & Veilleux, S. 2003, ApJS, 148, 383 
Norman, C. A., & Ikeuchi, S. 1989, ApJ, 345, 372 
Osterbrock, D. E., Fulbright, J. P., Martel, A. R., Keane, M. J., 

Trager, S. C. & Basri, G. 1996, PASP, 108, 277 
Rand, R. J., Kulkarni, S. R., & Hester, J. J. 1990, ApJ, 352, LI 
Rand, R. J. 1996, ApJ, 462, 712 
Rand, R. J. 1997, ApJ, 474, 129 
Rand, R. J. 2000, ApJ, 537, L13 



Reynolds, R. J. 1993, in AIP Conf. Proc. 278, Back to the Galaxy, 

ed. S. S. Holt & F. Verter (New York: AIP), 156 
Rossa, J., & Dettmar, R.-J. 2000, A&A, 359, 433 
Rossa, J., & Dettmar, R.-J. 2003a, A&A, 406, 493 
Rossa, J., & Dettmar, R.-J. 2003b, A&A, 406, 505 
Rownd, B. K. & Young, J. S. 1999, AJ, 118, 670 
Sancisi, R. & Allen, R. J. 1979, A&A, 74, 73 

Schaap, W. E., Sancisi, R., & Swaters, R. A. 2000, A&A, 356, L49 
Shapiro, P. R. & Field, G. B. 1976, ApJ, 205, 762 
Sofue, Y. & Rubin, V. 2001, ARA&A, 39, 137 

Swaters, R. A., Sancisi, R., & van der Hulst, J. M. 1997, ApJ, 491, 
140 

Takamiya, T. & Sofue, Y. 2002, ApJ, 576, L15 

Tullmann, R., Dettmar, R.-J., Soida, M., Urbanik, M., & Rossa, 

J. 2000, A&A, 364, L36 
Tullmann, R., Rosa, M. R., Elwert, T., Bomans, D. J., Ferguson, 

A. M. N., & Dettmar, R.-J. 2003, A&A, 412, 69 
Wakker, B. P. & van Woerden, H. 1997, ARA&A, 35, 217 
Wang, Q. D., Walterbos, R. A. M., Steakley, M. F., Norman, C. A., 

& Braun, R. 1995, ApJ, 439, 176 
Wolfire, M. G., McKee, C. F., Hollenbach, D., & Tielens, 

A. G. G. M. 1995, ApJ, 453, 673 
Wong, T. & Blitz, L. 2002, ApJ, 569, 157 



